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Abstract. A gradient-holonomic approach for the Lax type integrabiUty analysis of differential- 
discrete dynamical systems is described. The asymptotical solutions to the related Lax equation 
are studied, the related gradient identity subject to its relationship to a suitable Lax type spectral 
problem is analyzed in detail. The integrability of the discrete nonlinear Schrodinger, Ragnisco- 
Tu and Burgers-Riemann type dynamical systems is treated, in particular, their conservation 
laws, compatible Poissonian structures and discrete Lax type spectral problems are obtained 
within the gradient-holonomic approach. 



1. Preliminary notions and definitions 

Consider an infinite dimensional discrete manifold M C hC^', C™) for some integer m G Z+ and 
a general nonlinear dynamical system on it in the form 

(1.1) dw/dt^K[w], 

where w G M and K : M T{M) is a Frechet smooth nonlinear local functional on M and 
t G K is the evolution parameter. As examples of dynamical systems on a discrete manifold 
M C Z2(Z;C^) one can consider the well-known [5J[II] discrete nonlinear Schrodinger equation 

, . dUn/dt = i{Un+l - 2Un + M„-l) - ZW„U„(m„+1 + 

dVn/dt = -i{Vn+l - 2Vn + W„-l) + W„M„(u„+l + Vn-l) 

the so called Ragnisco-Tu [26 equation 

(1.3) j^""/f " " """^ ] Kr^M, 

^ ' dv^ / dt = -Vn-l + Unvl J "1 ' J' 

where we put w := {u, £ M, and the inviscid Riemann-Burgers equation ^20^ on a discrete 
manifold M C Z2(Z;R) : 

(1.4) dwn/dt = Wniwn+i - it;„_i)/2 := Kn[w] 

and its Riemann type [3TJ [55] generalizations, where w G M, having applications ^ in diverse 
physics investigations. 

For studying the integrability properties of difFerential-difFerence dynamical system (|1.1|) we will 
develop below a gradient-holonomic scheme before devised in [6l I15[ 1131 [7] for nonlinear dynamical 
systems defined on spatially one-dimensional functional manifolds and extended in [12j on the case 
of discrete manifolds. 

Denote by (•, ■) the standard bi-linear form on the space T*{M) x T{M) naturally induced by 
that existing in the Hilbert space Z2(Z;C™). Having denoted by 2?(Af) smooth functionals on M, 
for any functional 7 G V{M) one can define the gradient grad 7[w] G T*{M) as follows: 

(1.5) grad ^[w] := j''*[w] ■ 1, 

where the dash-sign ""' means the corresponding Frechet derivative and the star-sign " * " means 
the conjugation naturally related with the bracket on T*{M) x T{M). 
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Definition 1.1. A linear smooth operator -d : T*{M) x T{M) is called Poissonian on the manifold 
M, if the bi-linear bracket 

(1.6) {;■}, :=(grad(-),i?grad(.)) 

satisfies [U [21 111 UHl E] the Jacobi identity on the space of functionals 'D{M). 

This means, in particular, that bracket (jl.Sp satisfies the standard Jacobi identity on 'D{M). 

Definition 1.2. A linear smooth operator d : T*{M) x T{M) is called Notherian [HI IS [8] subject 
to the nonlinear dynamical system (jl.ip . if the following condition 

(1.7) Lk^ = ^' K - dK''* - K'd = 

holds identically on the manifold M, where we denoted by Lk the corresponding Lie-derivative 
[H El [g ig along the vector field K : M ^ T{M). 

Assume now that the mapping d : T*{M) x T{M) is invertible, that is there exists the inverse 
mapping '0~^ := fi : T*{M) x T{M), and called symplectic. It then follows easily from (|1.7p then 
the condition 

(1.8) Lxfi^nK + nK' + K'-*n^o 

hold identically on M. Having now assumed that the manifold M C Z2(Z;C^) is endowed with a 
smooth Poissonian structure i9 : T*{M) x T{M) one can define the Hamiltonian system 

(1.9) dw/dt -.^ --d giad H[w], 

corresponding to a Hamiltonian function H G V{M). As a simple corollary of definition (|1.9p one 
obtains that the dynamical system 

(1.10) - d grad H[w] K[w] 

satisfies the Notherian conditions (|1.7p . Keeping in mind the study of the integrability problem 
m 131 [T31 IHj subject to discrete dynamical system (|l.ip . we need to construct a priori given set of 
invariant with respect to it functions, called conservation laws, and commuting to each other with 
respect to the Poisson bracket (|1.5p . The following Lax criterion f3l [6l [13] proves to be very useful. 

Lemma 1.3. Any smooth solution ip G T*{M) to the Lax equation 

(1.11) Lk ip = dip/dt + K'^*ip^O, 
satisfying with respect to the bracket (•, •) the symmetry condition 

is related to the conservation law 

1 

(1.12) 7:= JdX{ip[wX],w). 



Proof. The expression (|1.12p easily obtains from the well-known Volterra homology equalities: 
11 1 1 

(1.13) 7 = J-^^^dX^ J dX{l,j'[wX] -w,) ^ JdX{j''*[wX]-l,w) ^ JdX{giad-f[wXlw) 

00 

and 

(1.14) (grad7H)' = (grad 7M)'^*, 

holding identically on M. Whence one ensues that there exists a function 7 G P(Af ), such that 

(1.15) LkJ — 0, grad 7[w] — ip[w] 

for any w G M. □ 
The Lax lemma naturally arises from the following generalized Nother type lemma. 
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Lemma 1.4. Let a smooth element ip e T*{M) satisfy the Mother condition 

(1.16) Lk^j ^ d^/dt + K'^*^P ^ gmd 

for some smooth functional £ I?(M). Then the following Hamiltonian representation 

(1.17) K = grad 
holds, where 

(1.18) i?:=V''-V''* 
and the Hamiltonian function 

(1.19) H^ = {^,K)-C^. 

It is easy to see that Lemma [T751 follows from Lemma II. 4[ if the conditions i/j — ip * and = 
are imposed on (|1.16l) . 

Assume now that equation ()1.16|) allows an additional not symmetric smooth solution (p G 
T*{M) : 

(1.20) LK(l) = d(j)/dt + K''*(t) = grad C^. 
This means that our system (jl.ip is bi-Hamiltonian: 

(1.21) - 1? grad H^ = K = -7] grad i?^, 
where, by definition, 

(1.22) r]:^(t>' -</>'*, H^ ^ (c^, K) ~ C^. 

Definition 1.5. One says that two Poissonian structures "d^r] : T*{M) — > T{M) on M are com- 
patible [ini m m [S], if for any A, ^ e R the linear combination X-d + ^rj : T*{M) T{M) will be 
also Poissonian on M. 

It is easy to derive that this condition is satisfied if, for instance, there exist the inverse operator 
'd-^ : T{M) T*{M) and the expression r]{d-^ri) : T*{M) T{M) is also Poissonian on M. 

Concerning the integrability problem posed for the infinite-dimensional dynamical system (|l.ip 
on the discrete manifold M it is, in general, necessary, but not enough [4l[6l[T3], to prove the exis- 
tence of an infinite hierarchy of commuting to each other with respect to the Poissonian structure 
(|1.5|) conservation laws. 

Since in the case of the Lax type integrability almost of ()1.1|) there exist compatible Poissonian 
structures and related hierarchies of conservation laws, we will further constrain our analysis by 
devising an integrability algorithm under a priori assumption that a given nonlinear dynamical 
system (jl.ip on the manifold M is Lax type integrable. This means that it possesses a related Lax 
type representation in the following, generally written form: 

(1.23) A/„ :=/„+! -/„[u;;A]/„, 

where / := {/„ £ C" : n G Z} C ?2(^; C^) for some integer r G Z+ and matrices ln[w; A] G EndC^ , 
n G Z, in (|1.23p are local matrix- valued functionals on M, depending on the " spectral" parameter 
A G C, invariant with respect to our dynamical system (II. ip . 

Taking into account that the Lax representation ()1.23p is 'local' with respect to the discrete 
variable n G Z, we will assume for convenience, that our manifold M := Af(Ar) C /cx)(Z;C'") is 
periodic with respect to the discrete index n G Zat, that is for any n G Zat := Ij/NZ, and A G C 

(1.24) I4w; X] ^ ln+N[w; X] 

for some integer G Z+. In this case the smooth functionals on M(jv) can be represented as 

(1.25) 7 := E ^"M 



for some local 'Frechet' smooth densities 7„ : Af(jv) C, n G Z 



N- 
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2. The integrability analysis: gradient-holonomic scheme 

Consider the representation (|1.23p and define its fundamental solution F,„.„(A) £ Aut{€!.^), 
m,n £ Zat, satisfying the equation 

(2.1) 

and the condition 

(2-2) F™,„(A)|„^„ = 1 

for all A G C and n e Z^r. Then the matrix function 

(2.3) 5'n(A) Fn^N^nW 

is called the monodromy matrix for the linear equation (jl.24l) and satisfies for all ?i G Zat the 
following Novikov-Lax type relationship: 

(2.4) Sn+l{X)ln=lnSniX)- 

It easy to compute that S'„(A) Y[ ln+k[w;X\ owing to the periodical condition ()1.24p . 

fc=6,jv-i 

Construct now the following generating functional: 

(2.5) 7(A) :=tr5„(A) 
and assume that there exists its asymptotical expansion 

(2.6) 7(A) ^ J2 

as A ^> cxD for some fixed jo G Z+ . Then, owing to the evident condition 

(2.7) -D„7(A) = 

for all n G Z^v, where we put, by definition, the 'discrete' derivative 

(2.8) A. := A - 1, 

we obtain that all functional 7j G X'(M(a')), j G '^+, are independent of the discrete index 
n and are simultaneously conservation laws for the dynamical system 

Assume additionally that the following natural condition holds: the gradient vector 

(2.9) (^„(A) := grad 7(A)„N = tr Z''*(^„(A)C'), 

which solves the Lax determining equation (jl.lip . satisfies, owing to ()2.4|) . for all A G C the next 
gradient relationship: 

(2.10) z{X)d ^(A) = ^(A), 

where z : C — )• C is some meromorphic mapping, § and 77 : T*(M(^]\[-^) — > T(A/(-Ar)) are compatible 
Poissonian on the manifold M(^jq) and Notherian operators with respect to the dynamical system 
()1.1|) . As a corollary of the above condition one follows easily that the generating functional 
7(A) G 2?(M(7v)) satisfies the commutation relationships 

(2.11) {7(A), 7(a^)}^ = = {7(A), 7(m)}, 

for all A, /i G C. Thereby, if to define on the manifold Af(jv) a- generating dynamical system 

(2.12) dw/dr := -i? grad 7(A)H 

as A — >■ 00, it follows easily from (|2.11l) that the hierarchy of functionals (|2.6p is its conservation 
law. 

Since the existence of an infinite hierarchy of commuting to each other conservation laws is 
characteristic concerning the Lax type integrability of the nonlinear dynamical system (jl.ip . this 
property can be effectively implemented into the scheme of our analysis. Namely, the following 
statement holds. 
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Proposition 2.1. The Lax equation allows the following asymptotical as A — > cxd periodical 

solution Lp{X) e T*[M(^^^) : 

n 

(2.13) (pn{X) - a„(A)exp[Lj(t; A)]]J(T,(A), 

3=0 

where for all n E Z 

(2.14) a„(A) : = (1, a(i)^„['i«; A], a(2)^„[w; A], a(„_i)_„[u;; A])"", 
a(fc).„(A) ^ ^ag„HA-^+^a,(A)^ ^4^)NA-^+^ 



k ~ l,m — 1 and Lo{t; •) : C — >■ C, i G M, is some dispersion function. Moreover the functional 
7(A) ■— lii(A~'^CT„[w; A]) S X'(M(^)) zs a generating function of conservation laws for the 

dynamical system 

Proof. Owing to Lemma (|1.3p and relationship (|2.9p functional (|2.5I) is a conservation law for our 
dynamical system (|l.ll) . Based now on expression (12.31) and equation (|1.23l) one arrives at the 
solution representation (I2.13|) to the Lax equation (jl.lip . Now, making use of the periodicity of 
the manifold M(jv) , we obtain from the period translation of (|2.13p that the functional 



(2.15) 7(A) := MA-V„[7«;A]) ^ ^ 7,A- 



generates an infinite hierarchy of conservation laws to (|l.ip . so finishing the proof. □ 



Thus, if we start the Lax type integrability analysis of a priori given nonlinear dynamical 
system (jl.ip . it is necessary, as the first step, to study the asymptotical solutions (|2.13p to the 
corresponding Lax equation (II. lip and construct a related hierarchy of conservation laws in the 
functional form (j2.15p . taking into account expansions (j2.14p . 

Remark 2.2. It is easy to observe that, owing to the arbitrariness of the period N G Z+ of 
the manifold M(^), all of the finite-sum expressions obtained above can be generalized to the 
corresponding infinite dimensional manifold M C Z2(Z;C'"), if the corresponding infinite series 
persist to be convergent. 

Since our dynamical system (jl.ip under the above conditions is a bi-Hamiltonian flow on the 
manifold M,-^-), as the next step of the related compatible Poissonian or symplectic structures, 
satisfying, respectively, either equality (II. 7p or equality (jl.Sp . Before doing this, we need to 
formulate the following lemma. 

Lemma 2.3. All functionals jj G P(Af(jv)), entering the expansion 12.15]) . are commuting to 
each other with respect to both Poissonian structures §,7] : r*(A/(jv)) ^ 2^(-^^(JV))j satisfying the 
gradient relationship 12.16]) . 

Proof. Based on representations (|2.13p and (|2.9p one obtains that there holds the asymptotical as 
A — > 00 relationship 

(2.16) ln7(A)^7(A). 

Since the generative function 7(A) e 2?(Af(7v)) satisfles the commutation relationships ()2.1ip . 
the same also holds, owing to ()2.16p . for the generating function 7(A) G X'(Af(7v)), flnishing the 
proof. □ 

Proceed now to constructing the related to dynamical system (jl.ip Poissonian structures rj : 
r*(Af(jv)) '^(-^^(A'))- Note here, that these Poissonian structures are Notherian also for the 
whole hierarchy of dynamical systems 

(2.17) dw/dtj :— ~d grad 7j[w], 

where tj G M, j G Z+, are the corresponding evolution parameters, and which, owing to ()2.1ip . 
commute to each other on the manifold Af(Ar). The latter makes it possible to apply Lemma 11.41 
to arbitrary one of the dynamical systems ()2.17p , if the related vector fields commuting with (jl.ip 
are supposed to be found before. 
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To solve analytically equation (|1.16p subject to an element Lp G T*(Af(^)) one can, in the case 
of a polynomial dynamical system (ll.ip . make use of the well known asymptotical small parameter 
method [1^1 IS]- If applying this approach, it is necessary to take into account the following 
expansions at zero - element = € M( ^r) with respect to the small parameter /i : 

d/dt = d/dto + ^d/dti + fi'^d/dt2 + 

(2.18) K[w] = /iif^i^H +Ai^^^if^^^N + 
K'[w] = + nK[[w]+ fi^K!^[w] + 

grad C[w] = grad + ^grad C^^'^ [w] + ^^grad C^'^'^ [w] + ... . 

Having solved the corresponding set of linear nonuniform functional equations 

VVrf^o + i^oV^"^ = grad/:("\ 

(2.19) dip^^^/dto + K'o*^'-^^ = grad £(1) - A-;;V^°\ 
dip^^^dto + K'*^^^^ = grad £(2) _ A'^Vd) - ^^^(o) 

and so on, by means the standard Fourier transform applied to the suitable iV-periodical functions, 
one can obtain the related Poissonian structure as the scries 

(2.20) = (^(")^' - </.(0)-'* + - (^(i)-'*) + ... 

and next to put in (|2.20p at the end /i = 1. 

Another direct way of obtaining an Poissonian operator d : T*{M(^j^)) — T(M(jv)) for (II. ip is to 
solve by means of the same asymptotical small parameter approach the Notherian equation (jl.7p , 
having preliminary reduced it to the following set of linear nonuniform equations: 

(2.21) :^{^0V^°^) = ifo(^o¥'^"^), 

-^(i?i^(°)) = if;^(^i^(")) + ^oif;^>^°^+^iW°\ 

A(^,<^(0)) = if;;(^2^("))-^(o)'i^i + ^„i^^V(°' + 
ato 

Based now on the analytical expressions for actions -dj : ip'-^^ 'djip'^^\ j G Z+, one can easily 
retrieve them in operator form from the expansion 

(2.22) i? = i9o + M^i+M^i?2 + ..., 

if to put at the end of calculations jj, — 1. Similarly one can also construct the second Poissonian 
operator rj : T*(M(jv)) T{M(^n)) for the nonlinear dynamical system (jl.ip . 

Resuming up all this analysis described above, we can formulate the following proposition. 

Proposition 2.4. Let a nonlinear dynamical system on a discrete manifold Af(-jv) allow both a 

nontrivial symmetric solution Lp £ r*(M(jv-)) to the Lax equation il.ll]) in the asymptotical as A — > 
oo form \2.13\) . generating an infinite hierarchy of nontrivial functionally independent conservation 
laws \2.15\) . and compatible nonsymmetric solutions "0 and (j) G T*{M(^pfj) to the Mother equations 
iL16\ and U.20\) . respectively. Then this dynamical system is a Lax type integrable bi-Hamiltonian 
flow on M(jv) with respect to two compatible Poissonian structures d,ri : r*(Af(jv)) ~^ 2^(-^^(Af))j 
whose adjoint Lax type representation 

(2.23) dA/dt = [A,if''*], 

where A := d~^ri is the so-called recursion operator, can be transformed, owing to the gradient 
relationship \2.1U\) . to the standard discrete Lax type form 

(2.24) dlnldt = [pn{l),ln] + (A.Pn(O)^n 

for some matrix Pn{l) S End<C^ describing the related to il.23\) temporal evolution 

(2.25) dfn/dt = pnil)fn 

for f elooi^;Cn. 
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Remark 2.5. Based on the property that ah Hamihonian flows (|2.17p commute to each other and to 
dynamical system (11.11) and using the fact that they possess the same Poissonian and compatible 
(t?, 77)-pair, the analytical algorithm described above can be equally applied to any other flow, 
commuting with ()1.1|) . 

Concerning the discrete linear Lax type problem (|1.23l) , it can be constructed by means of the 
gradient-holonomic algorithm, devised in [SJ [71 113| for studying the integrability of nonlinear dy- 
namical systems on functional manifolds. Namely, making use of the preliminary found analytical 
expressions for the related compatible Poissonian structures ??, 77 : r*(M(jv)) — > 7'(Af(Af)) on the 
manifold M(7v) and using the fact that the recursion operator A :— -d^^rj : r*(M(jv)) — ^ '^*i^{N)) 
satisfies the dual Lax type commutator equality (|2.23|) . one can retrieve the standard Lax type 
representation for it by means of suitably derived algebraic relationships. As a corollary of Propo- 
sition 12.41 one can also claim that the existence of a nontrivial asymptotical as A 00 solution 
to the Lax equation (II. lip can serve as an effective Lax type integrability criterion for a given 
nonlinear dynamical system (jl.ip on the manifold Mj^v). 

3. The Bogoyavlensky-Novikov finite-dimensional reduction 

Assume that our dynamical system (jl.ip on the periodic manifold Mj^vj is Lax type integrable 
and possesses two compatible Poissonian structures 'djrj : T*(M(7v)) T(M(7v)).Thus, we have 
the nonlinear finite-dimensional dynamical system 

(3.1) dwn/dt :— Kn[w] — — grad Hn[w] 

for indices n G Zin owing to its A^-periodicity. The finite dimensional dynamical system p.ip can 
be equivalently considered as that on the finite-dimensional space M(jv) — (C™)^ parameterized by 
any integer index n G Zjv, and whose Liouville integrability is of the next our analysis. To proceed 
with studying the flow (13. ip on the manifold M(jv) , we will make use of the Bogoyavlensky-Novikov 
[rei m reduction scheme H [H H H] . 

Let A(M(jv)) := ® ^•'{^{n)) be the standard finitely generated Grassmann algebra pi [51 [TB] 

of differential forms on the manifold Then the following differential complex 

(3.2) AO(M(^r)) A Ai(M(^r)) ^'...A A^M^N)) ^ A^+i(AV)) 4 

where d : A(M(^')) — > A(Af|-jv)) is the external differentiation, is finite and exact. Since the 
discrete 'derivative' D„ := A — 1 commutes with the differentiation d : A(M(^)) A(M(^)), 
[Dn,d] = for all n G Zjv, and for any element a G A°(M(7v)) 

(3.3) grad( ^ D„a„[u;]) = 0, 
one can formulate the following Gelfand-Dikiy type [TT] lemma. 

Lemma 3.1. Let C[w] £ A*'(M(jv)) be « Frechet smooth local Lagrangian functional on the manifold 
M(7v)- Then there exists a differential 1-form a^^-* € A^{M(^i^^), such that the equality 

(3.4) dCn[w] =< grad £„[w],dw„ > -|-L»„a£^' [w] 
holds for all n G Z^r. 

Proof. One can easily observe that 

fn r\ r 1 V^^ d£n[w] dCn[w] ■ 

(3.5) dCn[w] = 2^ < 75 — —,dwn+j >= 2^ < 75 — —,A'dWn >= 

N-1 r 1 fN-l \ 

= < ^A~^' > +Dn ^ < Pj,dWn+j > , 

where, by definition, 

(3.6) pk -^Y.^' 



rAREMA A. PRYKARPATSKY, NIKOLAI N. BOGOLUBOV (JR.), ANATOLIY K. PRYKARPATSKY, AND VALERIY H. SAMOYLENKO 



for fc = 0, — 1. Having denoted the expression 

A 

(3.7) grad£„H:=^A 
one obtains the result p.4p . where 



^-1 Ptr r 1 
^ aCn [w\ 



JV-1 

(3.8) Oin\M '■= X! <P]^dWn+j > 

3=0 

is the corresponding differential 1-form on the manifold M( jv) . Q 
Applying now the ^-differentiation to expression p.4p we obtain that 

(3.9) - Dnuji^^w] =< d grad Cn[w],Adwn > 
for any n G Z, where the 2-forni 

(3.10) J^^[w]:^da^^^[w] 

is nondegenerate on M(jv), if the Hessian df^C['w]/d'^Wn is also nondegenerate. 
Assume now that the submanifold 

(3.11) M(N) {grad ^f'H ^ 0; w e M(^n)} , 
where, by definition, the Lagrangian functional 

N-l 

(3.12) /:W :=_^^+^c,7„ 



with 7j G 'D{M), j — 0, N — 1, for some N G Z+, being the suitable nontrivial conservation laws 
for the dynamical system (jl.ip . constructed before, and cj G C, = 0, TV — 1, being some arbitrary 
but fixed constants. As a result of p. lip and (|3.9p we obtain that closed 2-form w^^' G A^(M(jv)) 
is invariant with respect to the index n G Z^r on the manifold M(^f^y Moreover, the submanifold 
p. lip is also invariant both with respect to the index n G Z^r and the evolution parameter i G M. 
Really, for any n G Z^r the Lie derivative 

(3.13) Lifgrad C^-^^ = (grad C^-^^YK + K''*{grad C'^^^) = 0, 

smce the functional C'-^^ G I?(M(jv)) is a sum of conservation laws for the dynamical system 
(|l.ip . whose gradients satisfies the Lax condition (|l.lip . Moreover, it is easy to see that if the 

Lie derivative Lk grad = 0,n G Z^v, at i = 0, then grad = for all < G M 

and n G Z^r. Thus, the Bogoyavlensky-Novikov reduction of the dynamical system p.ip upon the 
invariant submanifold M(jv) is defined completely invariantly. 

Now a question arises: how are related the dynamical system (jl.ll) . naturally constrained to 
live on the submanifold M(jv)j and the dynamical system (jl.ip . reduced on the finite dimensional 
submanifold M(jv) C M(jv). To analyze this reduction we will consider the following equality: 

(3.14) < grad Ci^^[w], K^[w] >= -Dnh^^'>[w], 

for some local functional h^^'^l'w] G A°(M), following from the conditions (13. 3p and (|l.lll) : 

(3.15) grad < grad df^[wlK„[w] (grad Clf^[w]y-*K„[w] + X;;*Hgrad 4P[w] = 

= (grad /:f )H)'if„N + X:;*Hgrad C^f^w] = L^^grad CiP[w] = 0, 

giving rise to (|3.14p . Since on the submanifold M(n) the gradient grad ci^\w] = for all 
n G Zjv, we obtain from p.l4p that the local functional G A'^(M(jv)) does not depend on 

index n G Z^r. 

The properties of the manifold M(n)i described above, make it possible to consider it as a 
symplectic manifold endowed with the symplectic structure w^^^ G A^(M(7v)), given by expressions 
(|3.8p and (|3.10p . From this point of view we can proceed to studying the integrability properties of 
the dynamical system (jl.ip reduced upon the invariant finite-dimensional manifold A/(jv) C M(^[^y 
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First, we observe that the vector field d/dt on M(jv) is canonically Hamiltonian [TJ [31 |J with 
respect to the symplectic structure uj^"^^ e A^(M(jv)) : 

(3.16) -i^u;(2)(w,p) =d/i(*'(u;,p), 

where h'^''-\w,p) :— h^*'' (w) , (jJ^^\w , p) ui^^^^lw] and {w,py G M(^-) are canonical variables in- 
duced on the manifold M(n) by the Liouville 1-form p.Sp . Really, from expression p.l4p one 
obtains that 

did_ < grad £^^^[w],dw„ >= —Dndhl^^[w], 
which, being supplemented with the identity (j3.9p in the form 



id_d< grad £-l^\w],dWn >= — U^Zjia;^^-* [w] 



entails the following: 



(3.17) -J < grad cf\w],dw„ >= -D„{dh(^'>[w] + i^ujI^H^]), 
dt 

Since grad = = Lk grad C[w\ identically on M(n), from p.l7p one obtains the 

resuh (|3l6l) . 

The same one can claim subject to any of Hamiltonian systems (|2.17p . commuting with (jl.ll) 
on the manifold M. Moreover, owing to the functional independence of invariants 7j 6 2?(M(jv)), 
j = 0, iV — 1, entering the Lagrangian functional p.l2p . we can construct the set of functionally 
independent functions /i*^^^ G 2?(M(jv)), J = 0,iV — 1, as follows: 

(3.18) < grad Clf^[wld„ grad 7,-„[u;] D^hi^'^lw], 



It is easy to check that these functions ft.'^-'-' G 'D{M(^jsi-^), j = 0, TV — 1, are invariant with respect to 
indices n G Zjy and commuting both to each other and to the Hamiltonian function /i^*^ G P(M(jv)) 
with respect to the symplectic structure cj^^^ G A^(M(7v))- Thus, if the dimension dimM(jv) = 2iV, 
the discrete dynamical system (jl.ip reduced upon the finite-dimensional submanifold M(^m) 

M(jv) will be Liouville integrable. If the set of conservation laws "fj G I?(M(jv)), j = 0,7V — 1, 
proves to be functionally dependent on M^^j^-^, the described scheme should be modified by means 
of using the Dirac reduction technique [Tl[8l[6] for regular finding the symplectic structure a;(^)[u'] G 
A^(M(-jv)) on invariant nonsingular submanifolds. 

4. Example 1: the differential-difference nonlinear Schrodinger dynamical 

system and its integrability 

The mentioned before discrete nonlinear Schrodinger dynamical system p.2p is defined on the 
periodic manifold Af(jv) C ^oo(Z; C^). Its Lax type integrability was stated in [51[TT1[T3] making use 
of the simplest discretization of the standard Zakharov-Shabat spectral problem for the well-known 
nonlinear Schrodinger equation. In this Section we will demonstrate the application of the gradient- 
holonomic integrability analysis, described above, to this discrete nonlinear Schrodinger dynamical 
system (II. 2|) . First, we will show the existence of an infinite hierarchy of functionally independent 
conservation laws, having solved the determining Lax equation (II. lip in the asymptotical form 
(|2.13p . The following lemma holds. 

Lemma 4.1. The functional expression 

(4.1) := (^^^^^^ exp[zt(2 - A - \-')]f{a,{\), 

is an asymptotical, as A — > oo, solution to the determining Lax equation 

(4.2) d<^„M -f /C^ Vn - 

for all n G Z^r with the operator K''* : T*(M(7v)) 2^*(-^^(-/v)) of the form: 



(4.3) K'^* = 



-l{A -I- A ^) ■ VnUn 

V -^Un{un+l+Un-^:) +z(A + A-i).«„«„ 
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where, by definition, 



are the corresponding asymptotical expansions. 

Proof. To prove this Lemma it is enough to find the corresponding coefficients of the asymptotical 
expansions (|4.4p . To do this we will consider the following two equations easily obtained from 
g^ . g31) and gH) : 



(4.5) 



and 



(l-i;„_iw„_i)/i„(l- E ai^'A-^-i)-i - 1 

sez+ 

-W„(m„+1 + -U„_i) + iVn{Vn+l + W„_l) X] 



( E a5f^A-^)i^-i^[-ln/i„ + ln(l- ^ ai^^ A-^-^)] + 4z( ^ aJf^A-^) + 



(4.6) 



iXh„+i{vn+iUn+i - l){ E a^^'j^A ")( ^ a^'j^A ") - E ""''A 

sez+ sez+ sez+ 



(?;„_iu„_i - 1)( E ai%^\-')K,{l - E ^7).'^A-^-i)-i - E air^A- 
sez+ sez+ sGZ+ 

+ J7 I] al^'U"" - to„(m„+i + U„-l) + TO„(W„+1 + W„-l) E Ok^A""' 



seZ4 



Having equated the coefficients of (|4.5I) at the same degrees of the parameter A £ C, we obtain 
step- by-step the functional expressions for /i„, tri*' and Un^ 



(4.7) 



hn = 
.(0) = 



,(2) ^ 



(1 -«„u„),al") =0,aW =/3, 

w„_iu„_i + i;„_iu„_2(w„ +u„-2) - iDl{liihn-i)t, 

i-^^n-l + - 1) + a^^2iW„-l(u„ + U«-2), 



.(1) n-i/ 



dhn/dt ^ iDn{Vn-lUn - VnUn-l), 

for all n £ Z, s G Z, or 

(4.8) cr^°^ = W„_iU„„i -I- w„_iu„_2(w„ -I- U„-2) - «£'J;(ln/i„-i)f, 

O-i^' = + (1 ^ W„-lU„-l)(l - t;„-2Mn-2) + /3w„_l(u„ + M„-2), 

and so on. Thus, having stated that the corresponding iterative equations are solvable for all 
s £ Z+, we can claim that expression (j4.1l) is a true asymptotical solution to the Lax equation 

□ 



Recalling now that the expression 

N-l N-1 

(4.9) 7(A) - E ln/z„ + E ln(l - ^ ai^^X-^-'] 



n=0 



n=0 
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as A — > oo is a generating function of conservation laws for the dynamical system (jl.2p , one finds 
that functionals 



N-l N-1 



(4.10) 70 = ^ln(l-«„M„),7o = -$]a("), 

n=0 n=0 



71 



72 



71=0 

N-1 



3 

n=0 



and so on, make up an infinite hierarchy of exact conservative quantities for the discrete nonlinear 
Schrodinger dynamical system (|1.2I) . 

Make here some remarks concerning the complete integrability of the discrete nonlinear 
Schrodinger dynamical system (|1.2p . First we can easily enough state, making use of the standard 
asymptotical small parameter approach [SI [T31 [7] , that the Nother equation (|1.7p on the manifold 
M(7v) possesses [TH [TT] the following exact Poissonian operator solution: 



(4.11) ,9, 



ihn 
-ihn 



n G Z^r, subject to which the dynamical the dynamical system (jl.2p is Hamiltonian: 
(4.12) d{u, vy /dt = grad H4u,v] 

on the periodic manifold M(^^^, where the Hamiltonian function 



(4.13) Hi} := ^ \nhl - ^ {vnUn-i + «„-!"„-) = 21n7o - |(7o + 7o)- 



n=0 n=0 



By means of similar, but more cumbersome calculations, one can find the second Poissonian oper- 
ator solution to the Nother equation (I1.7P in the following matrix form: 



v„D-^VnA -(1 + VnD-yu„)A-^ 



ri-1 N-l 

where the operation D^^{...) :— {■■■)k — S (■••)fc] is quasi-skew-symmetric with respect 

k=0 k=n 

to the usual bi-linear form on T*(Af(Ar)) x T(A/(jv)), satisfying the operator identity {D^^)* = 
-A-^D-^A, neZ. 

The Poissonian operators (j4.1ip and (|4.14[) are compatible, that makes it possible to construct 
by means of the algebraic gradient-holonomic algorithm the related Lax type representation for 
the dynamical system (|1.2p . The corresponding result is as follows: the discrete linear spectral 
problem 

(4.15) Afn^ln[u,V;X]fn, 

where / £ IqcC^', C^) and for n G Z 



(4.16) ln[u,v;X] 



v„ X^^ 



allows the linear Lax type isospectral evolution 

(4.17) dfn/dt^pnil)fn 

for some matrix p„(Z) S End C^, n € Z, equivalent to the Hamiltonian fiow 

(4.18) dfn/dt^ {H^JnU, 
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where {., is the corresponding to (I4.11[) Poissonian structure on the manifold M(^j^y The equiv- 
alence of (|4.11l) and (|4.18l) can be easily enough demonstrated, if to construct the corresponding to 
(|4.15|) monodromy matrix Sn{X),n G Z, for all A e C and to calculate the Hamiltonian evolution 

(4.19) iSn(X) = {H^, Sn{\)U = [Pn{l),Sn{X)], 

at 

giving rise to the same matrix p„(/) £ EndC^ , n G Z, as that entering equation (j4.17p . 

Thus, we have shown that the nonlinear discrete Schrodinger dynamical system (jl.2p is Lax 
type integrable bi-Hamiltonian flow on the manifold M(jv)- Since the solution if{X) G T*(M(^r)), 
constructed above, satisfies the gradient-like relationship 

(4.20) Ai9 (/7(A) = 77 ip{X) 

for all for A G C, we derive that the found above conservation laws are commuting to each other 
with respect to both Poisson brackets {., and {■,■}■&■ The latter gives rise to the classical 
Liouville integrability [H [TS] of the discrete nonlinear Schrodinger dynamical system (|1.2I) on 
the periodic manifold M(jv)- The detail analysis of the integrability procedure via the mentioned 
before Bogoyavlensky- Novikov reduction [TBI H] and the explicit construction of solutions to the 
dynamical system (|1.2[) are planned to be presented in a separate work. 



5. Example 2: the Ragnisco-Tu differential-difference dynamical system and its 

integrability 

Consider the Ragnisco-Tu differential-difference dynamical system (|1.3p , defined on the periodic 
manifold M(jv) C ^oo(^;C^), and construct first the corresponding asymptotical solution to the 
Lax equation The following lemma holds. 

Lemma 5.1. The functional expression 

(5.1) ^„:=("'f))exp(Ai)na,(A), 

\ / j=i 

is an asymptotical, as A — > 00, solution to the determining Lax equation il.ll]) for all n ^ lif^ 
with the operator K' * : T*(Af(7v)) ^ ^*(-^^(-/v)) of the form: 

jr cy\ T^l.* _ ( A -'^ — 2UnVn v'^ 

where, by definition, 

(5.3) a„(A) ~ A(l- ^ ai'^[u,v]X'n, 

a„(A) ~ ^ a(f)[^.,^;]A-^ 

and the following analytical expressions 

(5.4) = -2un^iVn-2-D-'{dc7l^^/dt + a^^Ua^/^/dt), 

= -da^^^dt - 2{Un-lVn-2vl - UnVnvl_i), 



ana so on, 



hold. 



Proof. It is easy to calculate that local a- and a-functionals on A/( jv) satisfy the following functional 
equations: 

/g A(l - (T„(A)) + D-^jjlnaniX) - <a„(A) + 2u„w„ = 0, 

^ ■ ' da„(A)/dt + Aa„(A) +a„(A)D-i Jlna„(A) - 2w„w„A"ia„_i(A)cr„(A)-i 0, 



which allow the asymptotical as A — ?> oo solutions in the form (|5.3p . Then, solving step- by-step 
the corresponding recurrent equations one finds the exact analytical expressions ()5.4p . Taking now 
into account that for each n G there exists such a local functional Pn{X) that the expression 
J^lnCT„(A) ~ DnPn{X) holds on M(^) identically, we derive that the functional expression (|5.ip 
solves the Lax equation (|l.lip . proving the lemma. □ 
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As a simple corollary of Lemma 15.11 we obtain that the expression 

N 

(5.6) 7(A) := E HI E '^i:^^^'-') ^ E 7.^^' 

n=i sGZ-i- iez+ 

is a generating functional for the infinite hierarchy of conservation laws jj € Z3(M(7v)), J G of 
the Ragnisco-Tu differential-difference dynamical system (|1.3I) . 

Show now that the Ragnisco-Tu differential-difference dynamical system (|1.3I) is a bi- 
Hamiltonian dynamical system on the functional manifold M( jy) . Really, based on Lemma 11.41 
we can find that the element i/; := ^{vn,—Uny G T*(M(jv)) satisfies the functional equation 



N-l 

2 2 

2 



(5.7) dV-M + = grad £, = E 

fe=0 

giving rise to the first Poissonian structure 

1 



(5.8) K := ^'^ - i,'^* = 



-1 



on the manifold M(jv) with respect to which the differential-difference dynamical system ()1.3|) is 
Hamiltonian: 

(5.9) ■^{un,v„y = -'dn grad H^^n[u,v], 
where the Hamiltonian function, owing to relationship ()1.22p . equals 

(5.10) He := (^, K) -C^)^Y^ iulvl/2 - u„?;„_i) = E 



2 

fc=0 k=0 



-ul + 2u„D-iAm„ a - 2u„D-iAi;„ 



(2) 
n • 

The same way one can find the second compatible with (|5.8p Poissonian operator 

(5.11) 77„ := 
for which 

(5.12) ""-''^ ^ ^'^'^ ^^'^'^ -ff»),n[w, w], 
where the Hamiltonian function 

AT ^ N 

(5.13) H,, := - E ^riVn = 2 ^ "^"+1 ■ 

fe=l fc=l 

Moreover, we claim that the hierarchy of conservation laws (|5.6p satisfies as A ^ cx) the gradient 
relationship 

(5.14) \-d grad 7(A) — ?] grad 7(A), 

entailing their commutation with respect to both Poissonian structures (|5.8I) and (|5.11l) . The 
latter allows us to argue that the Ragnisco-Tu differential-difference dynamical system ()1.3p is a 
completely integrable bi-Hamiltonian dynamical system on the manifold M(^) . 

Since the gradient relationship (|5.14p gives rise to the following 'adjoint' Lax type representation 

(5.15) dA/dt = [A,if''*], 

where, by definition, the expression A :~ d^^rj : T*{M[js!^) — >■ T*{M(^i\fj) is called a recursion 
operator. Based on the gradient relationship ()5.14p and expression (|2.9p we can obtain, within the 
gradient holonomic approach, that the Ragnisco-Tu differential-difference dynamical system (|1.3p 
is also Lax type integrable whose standard linear shift Lax type spectral problem equals 

(5.16) A/„ = /„[u,w; A]/„, l„[u,v;X]^( ^ + 
for all ri G Z, AG C, where {u, v) G Afj^v) and / G Zoo(Z; C'^). 
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6. Example 3: the generalized differential-difference Riemann-Burgers and 
riemann type dynamical systems and their integrability 

In [29l |31] [30] there was analyzed by means of simple enough differential-algebraic tools a gen- 
eralized (owing to D. Holm and M. Pavlov) Riemann type hydrodynamical hierarchy of equations 

(6.1) D't U^O, Dt := d/dt + uD^,D^ :^ d/dx, 

on a smooth functional manifold Ai C C°° (R; R) for any integer s G and proved their both 
bi-Hamiltonian structure and Lax type integrability. For s = 2 equation (j6.1l) possesses a Lax 
type representation, whose Lax Z-operator is given [30] by the expression 



[u, v; A] 



2A2 u^X 



where the vector-function {u,v)'^ E A4 C C°°(]R;M^) satisfies the equivalent to (|6.ip nonlinear 
dynamical system 

(6.3) Dtu = V, Dtv = 0, 

and A e C is an arbitrary invariant spectral parameter. For studying differential-difference versions 
of equations (|6.3p we will apply to the Lax Z-operator (|6.2p the Ablowitz-Ladik discretization 
scheme [5]. As a result of simple calculations we obtain the following new discrete Lax type 
spectral problem: 

(6.4) /„+i = L[u, V- A]/„, l^[u, v-\]=l^^- 2^^"^" ^ + Ai£u„ 

for n G Z, where function / G /oo(Z; C^), the vector (u, vY E M, if the resulting discrete dynamical 
system is considered on an A^-periodical discrete manifold M C (K^)^™ . 

To study the related with the spectral problem (|6.4I) nonlinear differential-difference dynamical 
systems, we will make use of a slightly generalized gradient-holonomic scheme. 

First, we formulate the following simple but useful lemma. 

Lemma 6.1. The following matrix 

.g5^ p (^.( e^^UA) -Anzi™(A)g(^!„(A)/(2A) \ 

I -2Xii^liX)/DmiiraiX) e~^|„(A) i ' 



where 



2(1) 



n(A) : = Y[[l-XDkUk + X{DkUkf/DkilkW], 



k—7i 
m—1 

(6.6) gi'L(A) : =l[[l + X{DkUk^DkUk{m 

k—n 

jointly with functional relationships 



(6.7) 



_|_ D„M„(A)-D„M„ i-i 

,D^Un{X)J \d^u^ + A-i[l + """"^t^f)""" ]-^ 
and 

(6.8) 2DnVn - {DnU,,f = 0, 

solves the associated with spectral problem {6.4-^ linear matrix equation 

(6.9) Fm+1,„(A) = 
under the initial condition 

(6.10) F„,,„(A)U=„ =I+0(1/A) 
as A oo for all m,n £ 
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As a corollary from relationships (j6.7l) one easily finds limits 

(6.f 1) lim DnUn{X) = DnUn, lim DnUn{X) = D„Un, 

A— j-oo A— >oo 

uniformly holding for all n G Z. 

Now, based on the matrix expressions (j6.5p and (I6.6p . one can construct the fundamental matrix 

(6.12) F„,„(A) F^,„(A)F„-i(A), 
solving the linear problem 

(6.13) -Fm+l,n(A) = lm[u, V, A]i^™,„(A) 

under the initial condition 

(6.14) F„,„(A)| 
for all n £ Z as A — > oo. 

Taking into account that the manifold M is N - periodic, one can construct the next important 
objects - the asymptotical as A — > oo monodromy matrix 

(6.15) S„{X) F„+jv,„(A) 

for any n £ Z. By construction, the monodromy matrix (|6.15p satisfies the following useful prop- 
erties: 

(6.16) S„+n{X) ■■= SniX), det5„(A) = l, 

holding for all n e Z and A ^ cxo. 

Keeping now in mind the importance of invariants and Poissonian structures related with the 
linear spectral problem (j6.4p . we proceed to studying its basic Lie-algebraic properties and con- 
nections with the so called vertex operator representation [3H |33] of the related whole hierarchy 
of integrable differential-difference dynamical systems on the manifold M. 

Namely, we will sketch below the Lie-algebraic aspects [Ml [331 [SSI of differential-difference 
dynamical systems, associated with our Lax- type linear difference spectral problem (16.41) . where 
one assumes that the matrix In := ln[u, v; X] e G„ := GL2(C) (8) C(A, A~^) for n £ Z^v := Z/iVZ 
as A — >■ oo. To describe the related Lax type integrable dynamical systems, define first the matrix 

product- group :— ® Gj and its action x M^^^ ^^cf'' ^^'^ phase space M^^' :— 

{In £ Gn ■ n £ Zjv}, given as 

(6.17) {5„ £ G„ : n £ Zn} x {?„ £ G„ : n £ Zn} = {gnlngn+i £ G„ : n £ Zat}. 

Subject to action ( I6.17P a functional 7 £ ^(Mq^') is invariant iff the following discrete relationship 

(6.18) gTad-f{ln)ln = /„+igrad7(/„+i) 
holds for all n £ Zjv- 

Assume further that the matrix group G^ is identified with its tangent spaces Ti{G^), I £ G^, 
locally isomorphic to the Lie algebra G''^\ where is the corresponding Lie algebra of the Lie 
group G^ , isomorphic, by definition, to the tangent space Te{G^) at the group unity e £ G^ . 
With any element I £ G^ there are associated, respectively, the left 77; : Q^^^ Ti{G^) and right 
Pi ■■ 5^^' Ti{G^) diflFerentials of the left and right shifts on the Lie group G^, and their adjoint 
mappings : T*{G^) G^'^^-* and r]* : T*{G^) gi'^)-*, where 

(prgrad7(0, X) = (grad7(0, ^0 = (l grad7(0, X) tr(Z grad7(0^), 

(6.19) (r/f grad7(/), X) = (grad7(/), IX) = igTadj{l)l, X) trigmd-fil)lX) 

for any X £ Q^^^^ and smooth functional 7 £ V{G^), tr : G^ C is a trace-operation on the 
group G^ : trA :— resx=oD J2 SpAj[u, v; X] for any A £ G^ . Owing to (|6.18p and (|6.19p we can 

define the set 

(6.20) {$„ = grad7(Z„)^n S := r;(G), n £ Zjv} 

belonging to the space ~ T*{G^) and satisfying the following invariance property: 

(6.21) $„+l = Adl^niX) = ln'^n{X)ln 
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for any n G Zjv. The relationship (j6.2ip ahows to define a function ip : C invariant with 

respect to the adjoint action 

(6.22) G„ X G„ 9 [g, 5„(A)) ^ ad<,5„(A) = gSn{\)9-^ G G„ 
for any n ^"Ln and such that 

(6.23) 7(/) = p[Sn{\)]. = grad^[5jv(A)]5A,(A), 
where, by definition, the expression 

N 

(6.24) Sn{X) = ^U[u,v;X] 

i=i 

coincides exactly with the proper monodromy matrix for the linear spectral problem (16.41) . Since, 
owing to (I6.2ip . the matrices — gradiy9[S'„(A)]S'„(A) £ 5*, n E Z^r, can be reconstructed from 
(|6.24p . we find |M1I3S] the following Poissonian flow on the matrices S'„(A) e G„,n G Zjv : 

(6.25) dSn{X)/dt = [7^(gradv7[5„(A)]5„(A)), 5„(A)] 

with respect to the invariant Casimir function ip G I{Gn) and the quadratic Poissonian structure 

(6.26) {71,72} [l. [grad7l(0,7^(^ grad72(0)] + Ml grad7i (0), grad72(/)]) 

for any functionals 71,72 G 27(G^), constructed by means of a skew-symmetric 7?.-structure TZ : 
giN),* _^ g{N)^ particular, the equality 

(6.27) [grad(p(5„),5„] =0 
holds for all n G Zjy. 

Having taken into account (I6.23|) . one can rewrite (j6.25p in the following equivalent form: 

(6.28) dS„/dt = [7^(grad7(^„)^„), S„], 

holding for all n G Zjy. The latter jointly with (|6.2ip makes it possible to retrieve [551 [3S] the 
related evolution of elements In G G„, n G Z^r : 

(6.29) dln/dt = P,i+l{l)ln - lnPn{l), 

Pn{l) ■■ = 7^(grad7(Z„);„) 
and following from the relationships 

(6.30) ^„(A) = Vn(05iv(A)V^-i(0, 

ipn{l) = Illj[u,v;X\. 

Subject to the linear spectral problem (|6.4p the solution / G /oo(Z, C^) satisfies the associated 
temporal evolution equation 

(6.31) dfn/dt = pn{l)fn 

for any n G Z. It is easy to check that the compatibility condition of linear equations (|6.4p and 
(j6.3ip is equivalent to the discrete Lax type representation (j6.29p . which, upon reducing it on 
the group manifold Mg, gives rise to the corresponding nonlinear Lax type integrable dynamical 
system on the discrete manifold M. It follows from the fact that all Casimir invariant functions, 
when reduced on the manifold Mg, are in involution [35JI3B] with respect to the Poisson bracket 

O- 

In the case when the skew-symmetric 7^-structure TZ — 1/2(P+ — P_), where P± : Gn Qn,± C 
Qn are the projectors on the A-positive and A- negative, respectively, degree subalgebras of the Lie 
algebra Qm the determining Lax type equation (|6.29l) generates the flows 

(6.32) Az„[u,z;;A] - {X^ + ^Sn+l{\)) + ln[u,V-\]-ln[u,V-X]{X^ + ^Sn{\)) + 
dtj 

for all j G Z_|_ , where Sn (A) , n G Z^v , are the corresponding asymptotical expansion of the suitably 
normalized monodromy matrix S'„(A) E Q_,n E Zjv, as A — > 00. 

The hierarchies of evolution equations ()6.32p can be rewritten as the following generating flows: 

(6.33) -^ln[u,v;\] = [Sn+i{^i)ln[u,v; A] - l„[u,v; \]S„{n)] 
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as A — )■ oo and |A//x| < 1, where, by definition, 

d 



.34) 



dt 



dti 



Proceed now to describing the analytical structure of the regularized matrices 5,1 (/i), n € Zjv, for 
arbitrary /i g C. To do this effectively we need to consider the corresponding to (I6.33P evolution 
equations for the monodromy matrix S'„(A) as A cxd, which make it possible to construct the 
related flows on functions d„ := DnUnW and a„ := Z?„M„(A),n € Za?, represented [321 [33] by 
means of the related vertex operators action X\ : — M^, where M := M^", 



(6.35) 

and 
(6.36) 



{expDx,cxp{-Dx)y , 
1 d 



E 



(jA^)dt/ 



u(io + l/A,ti + l/(2A2),...) 
u(to-l/A,ti-l/(2A2),...) 



as A — )> oo. 

Now, based on the flows 16.331 one can derive the corresponding evolution equation on the 
monodromy matrix 5„(A) 6 G-,n G Z^, with respect to the vector field (|6.34l) : 



(6.37) 



dt 



(m) 



A' 



as |A//i| < 1, A oo, which entails upon taking the limit /i — s- A the equation 
(6.38) 

where we put, by definition, 
(6.39) 



d 
It 



dti 



Now, taking into account Lemma 16. 1[ the related matrix expressions (|6.5p . (|6.12l) and (|6.15p 
and having analyzed the analytical structure of the resulting monodromy S'„(A) € Q-,n G Z^r, one 
can state by means of simple but slightly cumbersome calculations the following proposition. 



Proposition 6.2. The following differential relationships 



(6.40) 



dt 



dt 



2 Ask 



1- Wl- 



,(11), (11) 



_ \2 d 

- dX 



^ dX 



2As<") 



2As 



1- Wl 



jii) jii) 



jointly with equalities 



.41) 



D„Un{X) 



2A, 



= (11) 



DnU„ (A) 

hold as A — 7- cx) /or all n G Z^r. 



2Ask 



1- Jl 



jii) jii) 



As a corollary of (|6.1ip and differential relationships (|6.40p one obtains easily the following 
vertex operator representation (|6.36l) for the functions (|6.41l) : 



(6.42) 



DnUn{X) 
AiWn(A) 



AiUn(io + l/A,ti + l/(2A2),...) 
i?nU«(to-l/A,ii-l/(2A2),...) 



= Xx 
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holding as A — ?► oo for all n £ Z^r. Recalling additionally the algebraic relationships (|6.7I) 
and (|6.8|) . one finds from (|6.42l) the related infinite hierarchy of differential-difference dynamical 
systems on the manifold M : 

d ( DnUn \ _ ( Cj{Dlun)~^ 



for all j e Z+, where Cj € M, j G Z_|_, are some recurrently calculated constant coefficients. 

Thus, we see that the constructed hierarchy (|6.43p does not contain a "naive" discretization of 
the generalized Rieniann type system of equations ()6.3|) in spite of the fact that the difference Lax 
type linear spectral problem ()6.4|) regularly reduces to the corresponding linear spectral problem 
for dynamical system (|6.3p . 

Nonetheless, such a situation is not faced in the case of the inviscid discrete Riemann-Burgers 
dynamical system (ll.4p : 



(6.44) dWn/dt = WniWn+l - W„-l)/2 := Kn[w], 

defined on an A^-periodical discrete manifold AI C I-zC^n',^)- Following the gradient- holonomic 
scheme of studying the integrability of ()6.44|) . we first state the existence of an infinite hierarchy 
of conservation laws and the corresponding bi-Hamiltonian formulation. 
Consider, owing to Proposition , the determining equation (jl.lip 

(6.45) difjdt - [(A - A-i)w„/2 - {wn+i - w„-i)/2]</j„ = 
and its asymptotical solution (p e T*{M) in the form (|2.13p : 

(6.46) ^^ = \{a,[w-\l 

where n G Z and the local functionals <7j[w, A], j G possess as A — >• oo the following expansions 

(6.47) a,KA]^ J2 '^f'MA-^ 

sez+ 

Having solved recurrently the resulting functional equations 

(6.48) i:',7^(lno-„)t -t- (w„_i/o-„_i - w„+icr„)/2 + (w„+i - w„-i) = 0, 
one finds easily the infinite hierarchy (15.61) of conservations laws: 

(6.49) 70 = ^ (w„ -I- w„-i),7i = 0,72 = ^ [(w„ -I- w„-i)^ -t- w„(w„_i -I- w„+i)], ...,72^+1 = 

n=0 n=0 

for all j G Now, applying to the hierarchy of conservation laws the approach of Lemma 
11.41 one can find by means of slightly cumbersome and lengthy calculations the following pair 
lijij : T*{M) —> T{M) of compatible Poissonian operators on the manifold M : 

(6.50) ■■= w«(A - A~^)w„, ry„ := (w„-u;„+iA^ - w„w„_i A"^)(u;„ -t- w„_iA~^). 

In particular, there easily obtains the Hamiltonian representation of the Burgers-Riemann system 
(lOD : 

AT-l 

(6.51) dwn/dt = -i?„grad Hi), - ^ (u;„ -|- w„_i)/2. 

n=0 

Moreover, the first Poissonian structure of (j6.50p allows the continuous limit lim Wn ■= w(a;), if 

Ax^O 

n— f 00 

nAx := x G M, to the well known [20] correct form 

(6.52) i9:={wd + dw){w + d-^wd)/2. 

Making use the Poissonian pair ()6.50p one can retrieve within the gradient holonomic scheme a Lax 
type representation related to the inviscid discrete Riemann-Burgers dynamical system (|6.44p . 
whose /-operator is given by the following matrix expression: 

(6.53) ln[w; = ( I 
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for n G Z and A G C. Mention only that the higher flows, generated by the inviscid Burgers- Riemann 
dynamical system (|6.44p . have nothing to do with the generalized Riemann type hydrodynamic 
systems (|6.ip and their discrete approximations. Thereby, it is necessary to develop a different 
approach to constructing their proper integrable discrete Lax type representation, compatible with 
the related continuous limit. 

7. Conclusion 

The gradient-holonomic scheme of direct studying Lax type integrability of differential-difference 
nonlinear dynamical systems described in this work appears to be effective enough for applications 
in the one-dimensional case, similar to the case [6l [HI [131 (121 [H] of nonlinear dynamical sys- 
tems defined on spatially one-dimensional functional manifolds. This algorithm makes it possible to 
construct simply enough an infinite hierarchy of conservation laws as well as to calculate their com- 
patible co-symplectic structures. As it was also shown, the reduced via the Bogoyavlensky-Novikov 
approach integrable Hamiltonian dynamical systems on the corresponding invariant periodic sub- 
manifolds generate finite dimensional Liouville integrable Hamiltonian systems with respect to the 
canonical Gelfand-Dikiy type symplectic structures. As interesting examples the complete in- 
tegrability analysis of the nonlinear discrete Schrodinger, the Ragnisco-Tu and Burgers-Riemann 
type dynamical systems was presented. 

Subject to different not-direct approaches to studying the integrability of differential-difference 
dynamical systems on discrete manifolds it is worth to mention the works [191 UHl HH (HI [23l 
[28l (381 [24l l39l l40] based on the inverse spectral transform and related Lie-algebraic methods 
[25l [27l l38l [24l (391 l40] . where there are constructed a priori both conservation laws and the 
corresponding Lax type representations. Concerning these approaches, many of their important 
analytical properties were constructively absorbed by the gradient-holonomic scheme of this work 
and realized directly as an algorithm. 
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